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Abstract 

Several small-molecule CDK Inhibitors have been identified, but none have been approved for clinical use in the past few 
years. A new series of 4-[{3-hydroxybenzylamino)-methylene]-4H-isoquinoline-l,3-diones were reported as highly potent 
and selective CDK4 inhibitors. In order to find more potent CDK4 inhibitors, the interactions between these novel 
isoquinoline-l,3-diones and cyclin-dependent kinase 4 was explored via in silico methodologies such as 3D-QSAR and 
docking on eighty-one compounds displaying potent selective activities against cyclin-dependent kinase 4. Internal and 
external cross-validation techniques were investigated as well as region focusing, bootstraping and leave-group-out. A 
training set of 66 compounds gave the satisfactory CoMFA model (q^ = 0.695, = 0.947) and CoMSIA model (q^ = 0.641, 
1^ = 0.933). The remaining 1 5 compounds as a test set also gave good external predictive abilities with r^pred values of 0.875 
and 0.769 for CoMFA and CoMSIA, respectively. The 3D-QSAR models generated here predicted that all five parameters are 
important for activity toward CDK4. Surflex-dock results, coincident with CoMFA/CoMSIA contour maps, gave the path for 
binding mode exploration between the inhibitors and CDK4 protein. Based on the QSAR and docking models, twenty new 
potent molecules have been designed and predicted better than the most active compound 12 in the literatures. The QSAR, 
docking and interactions analysis expand the structure-activity relationships of constrained isoquinoline-l,3-diones and 
contribute towards the development of more active CDK4 subtype-selective inhibitors. 
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Introduction 

Cyclin-dependent kinases (CDKs), a family of serine /threonine 
protein kinases, play a central role in the growth, development, 
proliferation and death of eukaryotic cells [1,2]. There are more 
than 1 3 CDKs of which 1 2 different cyclin families have been 
identified up to now, and different CDK/ cyclin combinations are 
active during each phase of the cell cycle [3-.5]. Among these 
CDK/cyclin complexes, the D/CDK4 and E/CDK2 complexes 
have been greatly concerned [6]. In the Gl-S phase transition, the 
retinoblastoma susceptibility gene family of proteins (Rb) was 
phosphorylated by the D-type cychns (Dl, D2 or D3) in 
combination with CDK4 and cyclin E/CDK2 complexes. 
Phosphorylation of the Rb activated the E2F transcription factors 
and resulted in the transcription of genes required for DNA 
synthesis. This kind of function exerted by D/CDK4 and E/ 
CDK2 complexes is positively regulated by the mitogenic signaling 
pathways and negatively regulated by the cyclin-dependent kinase 
inhibitors (CKIs) [7-13]. Inhibition of cyclin-dependent kinases 
(CDKs) with small molecules has been suggested as a strategy for 
treatment of cancer, based on deregulation of CDKs commonly 



found in many types of human tumors. Selective CDK inhibitors 
such as CYC-202 [14] and BMS-387032 [15] targeting CDK2, 
and PD0332991 [16] targeting CDK4 have been under clinical 
evaluations. Recently, a series of novel isoquinoline- 1 , 3-(2H, 4//)- 
diones have been found to possess excellent selective inhibitory 
activity against the CDK4 [7,8]. 

The three-dimensional quantitative structure-activity relation- 
ship (3D-QSAR) models derived from the most widely used 
computational methods, CoMFA (comparative molecular field 
analysis) and CoMSIA (comparative molecular similarity indices 
analysis), could be used to guide rational synthesis of potent novel 
inhibitors and now aimed to elucidate the structural features 
required for CDK4 inhibitors. The best developed models have 
been duly validated by a systemic external validation, on the basis 
of which a set of twenty new potent molecules have been designed 
and predicted stronger activity than 12, the most active compound 
reported in the literatures[7,8]. 

Molecular docking techniques have been extensively used as an 
important tool in the discovery of new smaU-molecule drugs for 
targeting proteins. [17-21] Based on an idealized representation 
that a ligand makes every potential interaction with the binding 
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sites, docking uses an incremental construction algorithm to place 
flexible ligands into a fiiUy specified active site. Surflex-Dock is 
particularly successful at eliminating false positive results and 
therefore used to narrow down the screening pool significantly, 
while retaining a large number of active compounds. The binding 
interactions of the isoquinolinedione derivatives within the CDK4 
active sites were discussed. The step-wise description of method- 
ology used for 3D-QSAR analysis, molecular docking and 
designing of new CDK4 inhibitors is as shown in Figure 1. 

Materials and Methods 

Data Set 

All isoquinoline-l,3-(2//,4/;/)-dione derivatives and their biolog- 
ical activities were collected from literatures [7,8] (Figure 2). A 
total set of 8 1 molecules were randomly segregated into training 
and test sets comprising 66 and 15 molecules, respectively, based 
on the following rules: (i) Diversity of the molecules was very 
necessary to assess the statistical significance, (ii) To avoid any 
redundancy or bias in terms of structure features and activity 
range, the information of the selected compounds must be clear 
and concise, (iii) The most active and least active compounds 
should be included in the training set. The activities of the CDK4 
inhibitors were reported in IC50 and converted to pIC50 by 
taking Log (1/IC50) for the convenience. The activity range from 
4.6 to 8.6 log units of these compounds provided a broad and 
homogenous data set for 3D-QSAR study. In general, the spread 
of activity should cover at least 3 log units for a reliable 3D-QSAR 
model [22]. 

Molecular Modeling and Alignment 

The structures of the derivatives were sketched in SYBYL 8.1 
(Tripos, Inc., St. Louis, MO, USA) molecular modeling package 
and Gasteiger-Hiickel charges were assigned to the atoms of all the 
compounds. A good alignment is the most essential for the quality 
and the predictive ability of CoMFA and CoMSIA models [23], 
and common substructure, pharmacophore or docking overlaps 
can be available to align molecules [24,25]. The isoquinoline- 1 , 



3-(2H, 4H)-dione ring with structural rigidity was selected as the 
common substructure and the compound 12 with the strongest 
inhibitory activity as the template molecule (Figure 3). It can be 
seen that aU the compounds studied have similar active 
conformations. 

CoMFA and CoMSIA Setup 

Three-dimensional grid spacing was set at 2 A in the x, y, and z 
directions. The steric and electrostatic field energies were 
calculated using the Lennard-Jones and the Coulomb potentials 
[26]. For CoMFA method, a sp' hybridized carbon atom with a+1 
charge was identified as the probe atom to determine the 
magnitude of the steric and electrostatic field values, whose 
truncation was set at 30 kcal/mol [27-29]. 

The CoMSIA method, similar to CoMFA in terms of fields 
around the molecule, was based on the assumption that changes of 
ligands in binding affinities are associated with changes of 
molecular properties. Besides steric and electrostatic fields, three 
other different fields of hydrophobic, hydrogen bond donor and 
hydrogen bond acceptor are also calculated [30]. Moreover, a 
Gaussian function introduced in similarity indices makes it be 
calculated at all grid points, inside and outside different molecular 
surfaces. Equation used to calculate the similarity indices is as 
follows: 



probe,k 



Where, A is the similarity index at grid point q, summed over all 
atoms i of the molecule 7' under investigation. li^pmbe, k is the probe 
atom with radius 1 A, charge +1, hydrophobicity +1, hydrogen 
bond donating +1 and hydrogen bond accepting +1. W-^ is the 
actual value of the physicochemical property k of atom i. r-^^ is the 
mutual distance between the probe atom at grid point q and atom i 
of the test molecule, a is the attenuation factor whose optimal 
value is normally between 0.2 and 0.4, with a default value of 0.3 
[31]. 
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Figure 1. Step-wise description of methodology used for 3D-QSAR analysis, molecular docking and designing of new inhibitors for 
CDK4. 
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Figure 2. Chemical structures and IC50 values of the training and test set molecules. 

doi:10.1371/journal.pone.0093704.g002 



Partial Least Squares (PLS) Analysis 

For Partial Least Squares (PLS) analysis [32], the "leave-one- 
out" cross-validation metliod was first carried out to generate a 
cross-validated r [q ) value and the optimal number of compo- 
nents (ONG), based on the lowest standard error of prediction 
(SEP) which usually corresponds to the highest cross-validated 
squared coefficient [q ). To avoid over-fitting the models, a higher 



component was accepted only when the f diflFerences between 
two components was larger than 10% [24]. Region focusing was 
performed to maximize q' value by rotating the extracted principal 
components [33]. The q^ is a good indicator of the accuracy of 
actual predictions and a q value of 0.5 means halfway between no 
model and a perfect model [34]. Non-cross-validation was then 
executed to establish the final 3D-QSAR model after the optimal 
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predictive ability based on tlie statistics for the test set, it is not 
always sure tliat the model will perform well on a new set of data 
[35]. Therefore, an external test sets (''^,.^^) [36] was recommended 
for the estimation of predictive ability. Predictive values r^^.^^ were 
calculated as follows: 




Figure 3. Alignment of the compounds used in the training set. 

doi:1 0.1 371 /journal.pone.0093704.g003 

number of components was determined. The consequential final 
PLS models gave the conventional correlation coefficient (r^), 
standard errors of estimate (SEE) and F ratio between the 
variances of calculated and observed activities. The equation for 
SEE is given below 



SEE-- 



PRESS 



pred ' 



-\ -(PRESS /SD) 



Therein, SD is the sum of squared chfiFerences between the 
measured activities of the test set and the average measured 
activity of the training set. 

Several other statistics such as r^„^, r^o, R and k were calculated 
using the following equations, and 3D-QSAR models were 
considered acceptable only if they satisfy the following conditions: 
/,„>0.5, />0.6, [(/-r/)//] <0.1, 0.85<k<1.15 and r\^>0.5 
[36,37]. 



R = 



T,(yi-yo)(yi-yp) 
T,(yi-yo)J2(yi-ypf 



'-0 = 1 



E(j,->':°)^ 

2 



Jliy-yp) 
J2yyi 



Zyf 

yi''°=ky^ 



Where, n and c is the number of compounds and components, 
respectively, and PRESS is the sum of squared deviations between 
predicted and actual activity values for each molecule in the test 
set. 



Where, jc, and yp are the actual and predicted activities, 
respectively, yo and yp are the average values of the observed 
and predicted pICso values of the test set molecules, respectively. / 
is the non-cross-validated correlation coefficient from PLS process. 



External Validation 

is often a useful but not sufficient criterion for model 
validation. In many cases, a model with high and values were 
proved to be inaccurate. Even though a model may exhibit a good 



Molecular Docking 

Surflex-Dock in SYBYL 8. 1, using a patented search engine and 
an empirical scoring function to dock ligands into a protein's 
binding site [19], was apphed to study molecular docking in the 



Table 1. PLS results of CoMFA and CoMSIA models. 






CoMFA(before region focusing) 


CoMFA (after region focusing) 


CoMSIA 


PLS statistics 


LOO cross (j^/SEP 


0.544/0.727 


0.695/0.348 


0.641/0.496 


Group cross q^/SEP 


0.578/0.706 


0.711/0.334 


0.638/0.498 


Non-validated r'/SEE 


0.914/0.294 


0.947/0.185 


0.933/0.210 


F 


104.511 


1 39.423 


121.534 


rebootstrap 


0.916±0.023 


0.965 ±0.010 


0.928±0.015 


^bootstrap 


0.285 ±0.098 


0.152±0.071 


0.231 ±0.135 


Optimal components 


6 


6 


6 


Field distribution% 


Steric 


45.0 


47.9 


16.0 


Electrostatic 


55.0 


52.1 


25.2 


Hydrophobic 






18.9 


H-bond Donor 






10.1 


H-bond Acceptor 






29.8 





doi:l 0.1 371 /journa!.pone.0093704.t001 



PLOS ONE I www.plosone.org 



April 2014 I Volume 9 | Issue 4 | e93704 



QSAR, Docking of Isoquinolindinone CDK4 Inhibitors 



Table 2. The actual pICso, predicted pICso (Pred.) and their residuals (Res.) of the training and test set molecules. 





Compd. 


pICso 


CoMFA 




CoMSIA 




No. 


Actual 


Pred. 


Res. 


Pred. 


Res. 


1 


6.027 


5.905 


—0.122 


6.329 


0.302 
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6.099 


—0.009 


6.005 
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6.721 


6.801 
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6.579 
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6.556 


—0.1 14 


5 


6.638 
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—0.422 
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—0.1 75 


32 


5.791 
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6.132 
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0.223 
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5.901 
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—0.238 
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44 
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0.046 
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0.003 
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6.104 
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Table 2. Cont. 







Compd. 


pICso 


CoMFA 




CoMSIA 




No. 


Actual 


Pred. 


Res. 


Pred. 


Res. 


51 


6 886 


6 789 


_Q Qgy 


7 099 


0 213 


52 


7 569 


7 500 


— Q 059 
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0 044 


53 
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—0 T 50 
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0 1 80 
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55 
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4 834 


0 033 


4 748 


_0 053 


56 


6 495 


5 909 


_0 535 


6 001 


_0 494 


57 


6 000 


6 079 


0 079 


5 867 


_0 1 33 


58 


5 553 


5 604 


0 051 


5 713 


0 160 


59* 


5 569 


5 803 


0 234 


5 737 


0 168 


60 


5 523 


5 543 


0 020 


5 624 


0 101 


61 


4 790 


4 643 


_0 147 


4 996 


0 206 


62 


6 432 


6 507 


0 075 


6 670 


0 238 


63 


6 357 


6 503 


0 146 


6 609 


0 252 


64 


7 699 


7 499 


_0 200 


7 803 


0 104 


65 


8 000 


7 865 


_0 135 


7 787 


_0 213 


66 


7 523 


7 570 


0 047 


7 632 


0 109 


67 


7 699 


7 398 


_0 301 


8 004 


0 305 


68 


7 523 


7 723 


0 200 


7 422 


_0 1 01 


69 


7 699 


7 609 


_0 090 


7 865 


0 166 


70 


5 886 


6 003 


0 117 


5 599 


—0 287 


71 


6 252 


6 331 


0.079 


6 271 


0 019 


72 


6 569 


6 796 


0 227 


6 358 


_0 21 1 


73 


7 222 


7 301 


0 079 


6 933 


_0 289 


74 


4.917 


5.230 


0.313 


5.274 


0.357 


75* 


6.886 


6.558 


-0.328 


6.587 


-0.299 


76 


6.959 


6.960 


0.001 


6.789 


-0.170 


77 


6.086 


5.867 


-0.219 


6.372 


0.286 


78 


7.097 


7.109 


0.012 


6.903 


-0.194 


79 


6.602 


6.604 


0.002 


6.517 


-0.085 


80 


5.903 


5.913 


0.010 


6.043 


0.140 


81 


5.108 


5.140 


0.032 


5.006 


-0.102 



*Test set molecules. 

doi:1 0.1 371 /journal.pone.0093704.t002 
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Figure 4. Graph of actual versus predicted pICso of the training set and the test set using CoMFA (Left) and CoMSIA (Right). 
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Table 3. Results of the external validation for CoMFA and ColVlSIA models. 



correlation 



Parameters 


'^pred 


aslope 


intercept 


coefficient R 


/rslope 






CoMFA 


0.875 


1.021 


-0.025 


0.950 


0.986 


0.669 


-0.079 


CoMSIA 


0.769 


1.201 


-1.319 


0.943 


0.991 


0.631 


-0.100 



doi:1 0.1 371 /journal.pone.0093704.t003 

present paper. The crystal structure of CDK4 with hgand 1 GIH 
was retrieved from the RCSB Protein Data Bank [38]. A 
protomol, a computational representation of the receptor's 
binding cavity to which putative ligands are aligned, was generated 
automatically with a threshold parameter of 0.31 and a bloat 
parameter of 1 A, and composed of a collection of fragments or 
probe molecules such as CH4, N-H, and C = O that characterize 
steric effects in the binding pocket, hydrogen bond donor and 
acceptor groups, respectively. [39,40] AH the water molecules and 
sulfate salt in CDK4 IGIH (receptor) were deleted, and hydrogen 
atoms and Gasteiger charges were added [41,42]. AH of the eighty- 
one ligands were docked sequentially into the binding pocket of 
CDK4 using the parameters previously optimized. Surflex-Dock 
total scores are expressed in logio(A^) to represent binding 
affmities. The scores of 10 docked conformers of each isoquino- 
linedione derivatives were ranked in a molecular spreadsheet, and 
the highest total score was taken into consideration for ligand- 
receptor interactions. To visualize the binding mode between the 
protein and ligands, the MOLCAD (Molecular Computer Aided 
Design) program was applied to calculate and display the surfaces 
of channels and cavities, as well as the separating surface between 
protein subunits. MOLCAD program provides several types to 
create a molecular surface, in which the Robbin surfaces 
illustrating the secondary structure elements of the binding 
structure was applied to build the MOLCAD Robbin and 
Multi-Channel surfaces displayed with several potentials. Other 
parameters were established in default. 

Results and Discussion 



training and the test set molecules were listed in Table 2 and 
depicted graphically in Figure 4. For the CoMFA model after 
region focusing, the leave-one-out cross-validated value was 
0.695 (>0.5) and non-cross-validated value was 0.947 with an 
optimized component of 6, standard error estimate (SEE) of 0.185 
and i<" value of 139.423. Contributions of steric and electrostatic 
fields were 0.479 and 0.521, respectively. 

The CoMSIA model comprising all five descriptors gave a g'^ 
value of 0.641 and value of 0.933 with an optimized component 
of 6, standard error estimate (SEE) of 0.210 and F value of 
121.534. Contributions of steric, electrostatic, hydrophobic, 
hydrogen bond donor and acceptor fields were 0.160, 0.252, 
0.189, 0.101 and 0.298, correspondingly. 

External Validation Results 

The calculated results of the external validation were listed in 
Table 3. For CoMFA and CoMSIA models, the calculated y^prcd 
values were 0.875 and 0.769, with the slope (a) values of 1.021 and 
1.201 (close to 1), intercept {b) values of -0.025 and -0.039 (close 
to 0) and the correlation coefficient (R) values of 0.950 and 0.943 
(close to 1), respectively. The valid values of 0.669 and 0.631 
(>0.5) as well as high slope of regression lines through the origin 
{k) values of 0.986 and 0.991 (0.85<k< 1. 15) and the calculated 
[(/-r/)//] values of -0.079 and -0.100 (<0.1) were also 
obtained respectively. These external validation statistics revealed 
that both the CoMFA and CoMSIA models possessed high 
accommodating capacities and they would be reliable for 
predicting the pICso values of new derivatives. 



CoMFA and ColVlSIA Analysis 

The results taken from the PLS analysis were summarized in 
Table 1. The actual versus the predicted pIC-jo values for the 



CoMFA Contour Maps 

Figure 5. depicted the CoMFA steric and electrostatic contour 
plots for the most active compound 12. For the steric field, the 




(a) 




Figure 5. Std* coeff contour maps of CoMFA analysis in combination with compound 1 2. (a) Steric fields: green contours suggest regions 
where bulky groups increase activity, while yellow contours indicate regions where bulky groups decrease activity, and (b) Electrostatic fields: blue 
contours represent regions where electron-donating groups increase activity, while red contours highlight regions where electron-withdrawing 
groups increase activity. 
doi:1 0.1 371/journal.pone.0093704.g005 
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Figure 6. Std* coeff contour maps of CoMSIA in combination witKi compound 12. (a) Steric contour map. Green and yellow contours refer 
to sterically favored and unfavored regions, (b) Electrostatic contour map. Blue and red contours refer to regions where electron-donating and 
electron withdrawing groups are favored, (c) Hydrophobic contour map. White and yellow contours represent regions where hydrophilic and 
hydrophobic substituent are favored, (d) Hydrogen bond donor contour map. The cyan and purple contours indicate favorable and unfavorable 
hydrogen bond donor groups, (e) Hydrogen bond acceptor contour map. The magenta and red contours demonstrated favorable and unfavorable 
hydrogen bond acceptor groups. 
doi:1 0.1 371 /journal.pone.0093704.g006 
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The most potential derivatives 10-12 possessed corresponding 
electron-withdrawing aromatic groups (e.g. A^-pyrrolyl, 3-fiiryl and 
3-thienyl) at Ri, whUe the activities of compounds 31-32, 42 and 
44-46 bearing electron-donating substituents (e.g. piperidinyl, 
-NHAc, morpholinyl, jV,jV-dimethylformamido-) decreased signif- 
icantly. This also explained why compounds 1, 3-5, 7—16 and 19- 
24 with the corresponding methoxyl, hydroxyl, 2-furyl, 3-furyl, 3- 
pyridyl or 4-pyridyl group at R3 showed much more active than 
derivatives 25-60 and 74-75 with electron-donating substituent 
(1-piperidinylmethyl, substituted piperazinyl). A blue contour 
around the N-10 and C-11 positions emphasized the extreme 
importance of the electron-donating aminomethyl group. 



Figure 7. The binding mode between selected compound 12 
and the ATP pocket of CDK4 (PDB code: 1GIH). 

doi:1 0.1 371 /journal.pone.0093704.g007 

green contours represent regions of high steric tolerance (80% 
contribution) and the yellow contours (20% contribution) for 
unfavorable steric effect. The electrostatic field defined blue 
contours (80%) and red contours (20%) for electron-donating and - 
withdrawing substituents, respectively. 

In Figure 5a, one huge green contour around the Ri position 
revealed that bulky substituents at this site would benefit the 
activity, and two huge yellow contours near the Rj and R3 
positions suggested bulky groups at these sites unfavorable. This 
may explain the facts that derivatives 10-12 with relative bulkier 
groups (e.g. A^-pyrrolyl, 3-furyl and 3-thienyl) at Ri displayed the 
strongest activity, while derivatives 25-28, 31-32, 38-42, 44-45, 
47, 55, 57-60, 74 and 80 bearing a relative bulkier substituents 
(e.g. 4-methylpiperazinyl, -CH2-(1 -piperidinyl)) at R3 position 
showed a weak activity. Especially, derivatives 5, 7 and 8 with the 
corresponding substituent of chloro, bromo and iodo showed their 
activities in the following order of 5<7<8. 

One red contour near the Ri and one red around the R2 and R3 
in Figure 5b indicated an electron-withdrawing group favorable. 



CoMSIA Contour Maps 

The CoMSIA steric, electrostatic, hydrophobic, hydrogen bond 
donor and acceptor contours plots for the compound 12 were 
shown in Figure 6. The CoMSIA steric and electrostatic contour 
maps (Figs. 6a and 6b) were similar as the CoMFA steric and 
electrostatic contour maps (Figure 5). For hydrophobic field, white 
(20% contribution) and yellow (80% contribution) contours 
highlighted hydrophilic and hydrophobic properties, respectively. 
Hydrogen bond donor and acceptor fields take the cyan (80%) and 
purple (20%) contours for hydrogen bond donor, and the magenta 
(80%) and red (20%) contours as favor and unfavor for hydrogen 
bond acceptor, correspondingly. 

In Figure 6c, one yellow contour near the Ri region 
demonstrated that a hydrophobic substituent at this site would 
benefit the activity. Most of the active derivatives involved in 
present study possessed a hydrophobic group (e.g. 3-furyl, jV- 
pyrrolyl, 3-thienyl, bromo-, chloro-, iodo-, phenyl) at Ri, while 
those with only a hydrogen atom (e.g. 2, 3, 6 and 25-27) exhibited 
significantiy decreased potencies. Three pieces of white contour 
around the R2 position highlighted the hydrophilic properties of 
compounds 1—25 and 62-69. One purple contour near the Ri 
position in Figure 6d suggested a hydrogen bond donor group 
unfavorable. Therefore, the compounds 10-12, 49, 50 and 52-53 
with hydrogen bond acceptor oxygen or nitrogen atoms at Ri site 
exhibited better potencies. Both one red and two magenta 
contours located near the R] position in Figure 6e revealed that 
the hydrogen bond acceptor field was not very important for this 
site. The magenta contour near the R2 and R3 sites indicated 
hydrogen bond acceptor properties favorable. The hydroxyl 
groups at R2 could act as hydrogen bond acceptor at the same 
time. Therefore, the magenta contour confirmed the importance 
of the hydroxyl group at this region. Compounds 7-12 and 23 
bearing a hydrogen bond acceptor substituent (methoxyl, 4- 
pyridyl) at R3 showed the most activities. 



Bulky, electron-withdrawing 
hydrophobic and hydrogen 
bond acceptor groups favored 



Hydro; 
bindin 



gar bond 
g region 




Hydrogen bond 
binding region 



Minor, electron- withdrawing 
hydrogen bond acceptor groups 
favored 



Minor, electron- withdrawing 
hydropliilic, hydrogen bond 
donor and acceptor groups 
favored 



Figure 8. Structure-activity relationship revealed by 3D-QSAR and docking studies. 
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Molecule 



Substituent Predicted plCs, 





CoMFA 


CoMSIA 


OMe 


8.736 


8.368 


OEt 


9.407 


9.073 



12 3-furyl OH 

Dl \J COOH 

D2 1-^ COOH OEt 9.311 8.965 



D3 14 COOH OEt 8.909 8.798 



D4 \J/ COOH OEt 9.748 9.334 

NC^O, 

D5 \J COOH OEt 9.190 8.756 

D6 \J COOH OEt 9.330 8.895 



Cl3Cv.,0. 

D7 \J SOjNH, OMe 9.236 8.689 



D8 li SO2NH2 OMe 9.392 8.791 



D9 \J COCH, OEt 9.317 8.873 
V 

DIO \} COCH3 OEt 9.490 9.138 

Dll \J SO,H OMe 9.815 9.454 

Dll \J CONH2 NH2 9.820 9.247 

FjC^S 

D13 \J SO3H OMe 9.386 9.098 

D14 \J CONH2 NH2 9.456 8.982 



D15 SO3H OMe 8.855 8.519 



D16 "^^^ COCH3 OMe 9.235 8.780 

D17 \J COCH, OMe 9.237 8.869 

D18 \J SO2NH2 OEt 9.374 8.787 

D19 \J SO2NH2 OEt 9.011 8.753 

D20 \J SO3H OEt 8.995 8.767 



Figure 9. Structures and predicted pICSO values of newly designed molecules. 
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PLCS ONE I www.plosone.org 



10 



April 2014 I Volume 9 | Issue 4 | e93704 



QSAR, Docking of Isoquinolindinone CDK4 Inhibitors 



Molecular Docking Analysis 

Figure 7. illustrated the binding modes between compound 12 
and the ATP pocket. The carbonyl group at C-3 position acted as 
a hydrogen bond a('ccj)tor by forming a H-bond with tiie -XII 
group of Leu83 residue; the imino group at N-10 position served 
as a hydrogen bond donor and formed H-bond with the carbonyl 
group of Glnl31 residue; the hydroxyl group at site acted as 
both hydrogen bond donor and acceptor and formed two H-bonds 
with the carbonyl group of Asp 145 and the -NH group of Asnl32 
residues, respectively. The observations taken from Figure 7 
satisfactorily matched the corresponding CoMSIA hydrogen bond 
donor and acceptor contour maps. 

Summary of Structure-activity Relationship 

The structure-activity relationship revealed by 3D-QSAR and 
docking studies was illustrated in Figure 8. In short, the bulky, 
electron-withdrawing, hydrophobic and hydrogen bond acceptor 
groups at Ri position are favorable; the minor, electron- 
withdrawing, hydrophUic, hydrogen bond donor and acceptor 
groups at R2 position may benefit the potency; the minor, 
electron-withdrawing and hydrogen bond acceptor substituent at 
R3 position would increase the activity. The carbonyl group at C-3 
and the imino group at N-10 site were essential for binding to the 
ATP pocket of CDK4. 

Designing of Potent Derivatives 

Based on the structure-activity relationship revealed by the 
present study, twent) novel isoquinoline- 1, 3-(2//,4//)-dione 
derivatives were designed. These molecules were aligned to the 
database and their activities were predicted better than compound 
12 by the best CoMFA and CoMSIA models established 
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previously, especially D4, Dll and D12 showed 10 folds more 
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work. 
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Author Contributions 

Conceived and designed the experiments: PSJZ JMW GX. Performed the 
experiments: JZ HK JMW YG. Analyzed the data: PS JZ YG MY. 
Gontributed reagents/materials/analysis tools: PS JG GX. Wrote the 
paper: PSJZ JMW. 



12. Lu H, Schulze-Gahmen U (2006) Toward understanding the structural basis of 
cychn-dependent kinase 6 specific inhibition. J Med Chem 49: 3826—3831. 

13. Horiuchi T, ChibaJ, Uoto K, Soga T (2009) Discovery of novel thieno[2,3- 
d]pyrimidin-4-yl hydrazone-based inhibitors of Cyclin D1-CDK4: synthesis, 
biological evaluation, and structure-activity relationships. Bioorg Med Chem 
Lett 19: 30,5-308. 

14. McClue SJ, Blake D, Clarke R, Cowan A, Cummings L, et al. (2002) In viu-o 
and in vivo antitumor properties of the cyclin dependent kinase inhibitor 
CYC202 (R-roscovitinc). IntJ Cancer 102: 463^68. 

15. Misra RN, Xiao HY, Kim KS, Lu S, Han WC, et al. (2004) 
N-(cycloalkylamino)acyl-2-aminothia- zole inhibitors of cyclin-dependent 
kinase 2. N-[5-[[[5-(l,l-dimethylethyl)-2-oxazolyl]methyl]thio]- 2-thiazolyl]-4- 
piperidinecarboxamide (BMS-387032), a highly efficacious and selective 
antitumor agent. J Med Chem 47: 1719-1728. 

16. Toogood PL, Harvey PJ, Repine JT, Sheehan DJ, VanderWel SN, et al. (2005) 
Discoverv of a potent and selective inhibitor of cyclin-dependent kinase 4/6. 
J .Med Chem 48: 2388-2406. 

17. Wang J, KoUman PA, Kuntz ID (1999) Flexible ligand docking: a multistep 
strategy approach. Proteins 36: 1—19. 

18. Spitzer CM, Wellenzohn B, Laggner C, Langer T, Liedl KR (2007) DNA minor 
groove pharmacophores describing sequence specific properties. J Chem Inf 
Model 47: 1580-1589. 

19. Jain AN (2003) Surflex: fitlly automatic flexible molecular docking using a 
molecular similarity-based search engine. J Med Chem 46: 499-51 1. 

20. Jain AN (2007) Surflex-Dock 2.1: robust performance from l^and energetic 
modeling, ring flexibility, and knowledge-based search. J Comput Aided Mol 
Des 21: 281-306. 

21. Amburc PS, (iangAval RP, Sangamwar Al (2012) 3D-C)SAR and molecular 
docking analysis of biphenyl amide derivatives as p38 alpha mitogen-activated 
protein kinase inhil:)itors. .\lolec Div 16: 377-388. 

22. Cramer RD 3rd, Patterson DE, Bunce JD (1989) Recent advances in 
comparative molecular field analysis (CoMFA). Prog Clin Biol Res 291: 161- 
165. 

23. Cho SJ, Tropsha A (1995) Cross-validated R2-guided region selection for 
comparative molecular field analysis: a simple method to achieve consistent 
results. J Med Chem 38: 1060-1066. 

24. ZhengJ, Xiao G, Guo J, Zheng Y, Gao H, et al. (201 1) Exploring QSARs for 5- 
hpoxygenase (5-LO) inhibitory activity of 2-substituted 5-hydroxyindole-3- 
carboxylates by CoMFA and CoMSL\. Chem Biol Drug Des 78: 314-321. 



PLOS ONE I www.plosone.org 



11 



April 2014 | Volume 9 | Issue 4 | e93704 



QSAR, Docking of Isoquinolindinone CDK4 Inhibitors 



25. Afantitis A, Mrlagraki G, Sarimvcis H, Koutrntis PA, Markopoulos J, ct al. 
(2006) A novel simple QSAR model for the predietion of anti-HIV activity using 
multiple linear regression analysis. Molcc Div 10: 405^14, 

26. Murumkar PR, Giridhar R, Yadav MR (2008) 3D -quantitative structure-activity 
relationship studies on benzothiadiazepine hydroxamates as inhibitors of tumor 
necrosis factor-alpha converting enzyme. Chem Biol Drug Des 71: 363—373. 

27. Lan P. Huans? ZJ, SunJR, Chen WM (2010) 3D-QSAR and Molecular Docking 
Studies fin Fused Pyrazoles as p38alpha Mitogen-Activated Protein Kinase 
Inhibitors. IntJ Mol'Sci 11: 3357-3374. 

28. Hu R, Barbault F, Delamar M, Zhang R (2009) Receptor- and hgand-based 3D- 
QSAR study for a series of non-nuclcosidc HlV-1 reverse transcriptase 
inhibitors. Bioorg Med Chem 17: 2400-2409. 

29. Reddy RN, Mutyala RR, Aparoy P, Reddanna P, Rcddy MR (20 1 0) An analysis 
of hydrophobic interactions of thymidylate synthase with methotrexate: free 
energy calculations in\'olving mutant and native structures bound to metho- 
trexate. J Mol Model 16: 203-209. 

30. Klebe Q, Abraham U (1999) Comparati\'e molecular similarity index analysis 
(GoMSIA) to study hydrogen-bonding properties and to score combinatorial 
Ubraries. J Comput Aided Mol Des 13: 1—10. 

31. Sivan SK, Manga V (2010) Molecular docking and 3D-QSAR studies on 
triazohnone and pyridazinone, non-nucleoside inhibitor of HIV- 1 reverse 
transcriptase. J Mol Model 16: 1169-1178. 

32. Zhang N, Zhong R (2010) Docking and 3D-QSAR studies of 7-hydroxycou- 
marin derivatives as CK2 inhibitors. Eur J Med Chem 45: 292-297. 



33. Lindgren ¥, Geladi P, Rannar S, Wold S (1994) Interactive variable selection 
(IVS) for PLS. Part 1: Theory and algorithms. J Chemometrics 8: 349—363. 

34. Clark M, Cramer RD 3rd (1993) The probability of chance correlation using 
partial least squares (PLS). QSAR 12: 137-145. 

35. Roy K, Paul S (2010) Docking and 3D-QSAR studies of acetohydroxy acid 
synthase inhibitor sulfonylurea derivatives. J Mol Model 16: 951—964. 

36. Golbraikh A, Tropsha A (2002) Beware of q2! J Mol Graph Model 20: 269-276. 

37. Roy PP, Roy K (2008) On some aspects of variable selection for partial least 
squares regression models. Qsar Comb Sei 27: 302-313. 

38. Ikuta M/Kamaia K, Pukasawa K, Honma 'I\ Maehida 'I\ et al. (2001) 
CrystaUographic approach to idcntilication of cyclin-dcpcndcnt kinase 4 
(CDK4)-specific inhibitors by using CDK4 mimic CDK2 protein. J Biol Chem 

276: 27548-27554. 

39. Ruppert J, Welch W, Jain AN (1997) Automatic identification and represen- 
tation of protein binding sites for molecular docking. Protein Sei 6: 524—533. 

40. Holt PA, Chaires JB, Trent JO (2008) Molecular docking of intercalators and 
groove-binders to nucleic acids using Autodock and Surflex, J Chem Inf Model 

48: 1602-1615. 

41 . Muthas D, Sabnis YA, Lundborg M, Karlen A (2008) Is it possible to increase hit 
rates in structure-based virtual screening by pharmacophore filtering? An 
investigation of the advantages and pitfEiUs of post-filtering. J Mol Graph Model 
26: 1237-1251. 

42. Clark RD (2008) A ligand's-eye view of protein binding. J Comput Aid Mol Des 
22: 507-521. 



PLOS ONE I www.plosone.org 



12 



April 2014 | Volume 9 | Issue 4 | e93704 



